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J>v Abstract 



This paper proposes a novel, highly effective spectrum sensing algorithm for cognitive radio and 
whitespace applications. The proposed spectral covariance sensing (SCS) algorithm exploits the different 
statistical correlations of the received signal and noise in the frequency domain. Test statistics are 
computed from the covariance matrix of a partial spectrogram and compared with a decision threshold 
to determine whether a primary signal or arbitrary type is present or not. This detector is analyzed 
theoretically and verified through realistic open-source simulations using actual digital television signals 

m 

■ captured in the US. Compared to the state of the art in the literature, SCS improves sensitivity by 3 

qq | dB for the same dwell time, which is a very significant improvement for this application. Further, it is 



shown that SCS is highly robust to noise uncertainty, whereas many other spectrum sensors are not. 
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I. INTRODUCTION 

The scarcity of prime (frequencies below 2 or 3 GHz) spectrum is a decades-old problem in 
wireless communication, and will be for the indefinite future. Recognizing that a static allocation 
of spectrum over time and space is highly suboptimal - for example, often less than 10% efficient 
ID— 01 - there has been a flurry of interest in finding ways to adaptively allocate spectrum. 
Cognitive radio (CR) is a promising approach, whereby transmitter-receiver pairs find unused 
spectrum (white space) and use it for communication flU-O. The idea of spectrum reuse has 
received regulatory support in the form of the FCC white space ruling authorizing cautious reuse 
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of underutilized spectrum in the licensed TV bands 0. The IEEE 802.22 standard [H, lUl is 
attempting to formalize a solution that will meet FCC approval and allow communication in 
these bands in the near future. 

Reliable spectrum sensing - whereby devices determine whether another (e.g. licensed) signal 
is present - is fundamental to the success of cognitive radio. Primary (licensed) users have priority 
to use the channel and the secondary users can only use the resource when it is not occupied. 
Determining when the channel is occupied is quite difficult, because a potential transmitter must 
guarantee with high probability that no one in its transmission range could possibly be listening 
to a primary/licensed signal. This puts quite strict requirements on the detection capability^, and 
the secondary (i.e. cognitive) system is required to robustly sense the presence of primary signals 
at very low SNRs. For example, an IEEE 802.22 system must be capable of detecting digital 
TV signals at an SNR of -21 dB: when the noise is over one hundred times stronger than the 
actual signal. 

A. Spectrum Sensing: Related Work 

Accordingly, spectrum sensing research has been active, resulting in numerous proposed 
sensing algorithms, which are well summarized in ifTOl and 0T|. Studies on spectrum sensing can 
be classified in two categories: blind detection and feature detection. Blind detection is universal 
since it does not rely on prior information such as signal characteristics, channel and noise 
power, but the performance is generally poor. Feature detection senses specific characteristics 
of a known signal, and typically shows better performance than blind detection at the expense 
of not being applicable to all possible primary signals. The simplest (blind) sensing algorithm 
is an energy detector but it suffers from severe degradation under uncertain noise power [fT2|- 
lfl4l . A more robust blind sensing method called covariance absolute value (CAV) detection |[6l 
exploits the uncorrelated nature of the noise, whereas the primary signal is correlated. While the 
detection performance is comparable to the energy detector for a given noise power, it is much 
more robust to uncertain noise power than the energy detector. 

There have been several studies on the detection of digital TV signal features, which is one 
of the primary services to be protected in the IEEE 802.22 systems, and hence a logical starting 
point for feature detection sensing algorithms. A Field Sync Correlation Detector (FSCD) and a 

'The terms detection and sensing are used interchangeably throughout this paper. 
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Segment Sync Autocorrelation Detector (SSAD) are proposed in 11151 . which exploit the repetitive 
nature of synchronization sequences but their performance is similar to blind detectors. Cordeiro 
et al. proposed two FFT based pilot detection methods, one sensing the pilot energy and the 
other sensing the pilot location [2]. Pilot location detection is robust to the noise uncertainty 
while pilot energy detection is not. Under a certain noise power, both the methods achieve fairly 
good performance, the best in the literature to date as far as the present authors can judge. 

The idea of signal detection using cyclic spectral density (CSD) or spectral correlation function 
(SCF) is introduced in |[T6l . The cyclostationary nature of the modulated signal generates unique 
cyclic frequencies that are known to the detector a priori. FFT based pilot location detector Q is 
similar to the CSD detector in the sense that they both detect frequencies that contain relatively 
strong power. 

Tandra and Sahai developed the concept of a useful theoretical limit called the SNR wa n on 
reliable sensing under uncertain noise power ifTTl . lfT8l . They also proposed a matched filter 
with run-time noise calibration [14] to combat the time-varying nature of the wireless channel. 
Although the matched filter is an optimal detection algorithm, it requires perfect synchronization 
and parameter estimation, which is generally impractical. 

B. Contributions and Organization 

In this paper, we propose a novel spectral covariance sensing (SCS) algorithm that exploits 
different statistical correlations of the signal and noise in the frequency domain. The SCS detects 
spectral features for maximum sensitivity, but also is applicable for non-flat spectrum signals. We 
analyze its performance theoretically and then verify those results through extensive simulations 
using actual DTV signals captured in a real environment. Our simulations are publicly available, 
so that our performance claims can be externally verified |fT9l , which seems particularly important 
for this application 11201 . 

Rigorous comparisons with the FFT based pilot energy detector [2j and the CAV detector |[6l 
show that SCS achieves 3 dB better sensitivity than the FFT based pilot detection method (and 
an even larger gain vs. CAV) with the same sensing time or equivalently, achieves the same 
sensitivity in 20% of the sensing time. In the extremely low SNR regime, a 3 dB gain is a 
very significant improvement. It is also shown that the sensitivity of the CAV can be improved 
by a lot when combined with parts of the proposed algorithm. Furthermore, we show that the 
technique is unusually robust to noise uncertainty. The SCS achieves the best sensitivity of any 
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spectrum sensing approach thus far proposed in the literature, to the best of our knowledge. 

The rest of this paper is organized as follows. The SCS algorithm is presented in Section 
HO Section [III] analyzes detection performance of the SCS, especially in the low SNR regime. 
The theoretical results are verified through simulations using actual data and requirements of the 
IEEE 802.22 standards in Section [IV] Detection performance with various parameter selection 
is explored and comparisons with the previous methods are also presented. Section IVl concludes 
the paper. 

II. Spectral Covariance Sensing 

Assume that there is either one or no primary transmitter to detect, and that the secondary 
node can be located inside or outside of the primary cell boundary. The detection problem can 
then be formulated as a binary decision under two hypotheses as in fl6], [fl4l : 

H : z(n) = w(n), (1) 
Hi : z(n) = s(n)+w(n), (2) 

where z(n) is the equivalent received signal at baseband, s(n) is the signal component of received 
samples and w(n) is the noise component. The signal component s(n) may have a DC term 
(sd(n) — s d) an d an AC term (s a (n)). Therefore the signal samples are assumed to be independent 
and identically-distributed (i.i.d) with mean E[s(n)] = and variance a 2 s = E[s a (n) 2 ]. The noise 
samples w(n) are zero-mean i.i.d white Gaussian noise with variance a, 2 . 

For simplicity, we assume that the primary signal has a flat spectrum with a pilot tone at 
low frequency, consistent with the ATSC signal BU, lfT31 . However, the proposed algorithm can 
also be applied to signals with non-flat spectrum. The SCS algorithm exploits correlation of 
spectral feature in the band of interest where the signal specific spectrum that distinguishes 
itself from other signal or noise is located. Once the spectrum of the received signal is obtained 
by periodogram estimation EE], its correlation is computed using the sample covariance matrix. 
The primary signal typically has unique non-flat spectrum, so it is highly correlated, whereas 
the noise spectrum is flat and is almost entirely uncorrected. To find the relative correlation, the 
autocovariance of the spectrum is compared with the total covariance. The signal is detected if 
the spectral correlation is higher than the predefined threshold. The proposed SCS algorithm is 
described as the following steps. 
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1) Downconvert the received signal x(t) to a baseband complex signal y(t) = x{t)e~^^ ct to 
place the pilot tone at or near DC, where f c is the pilot frequency. 

2) Low pass filter and downsample y(t) by appropriate sampling rate F s to form z(n). 

3) Compute ^(n)'s spectrogram by calculating the squared magnitude of its short-time Fourier 
transform (STFT) as 

1 



Z T {k) 



N 



N-l 
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(3) 



where N is the number of FFT points, r e {0, 1, • • • , Nd — 1} is the index of the sensing 
window, Nd is the total number of sensing windows and k E {—N/2, ■ ■ ■ , 0, • • • , N/2 — 1} 
is the frequency index. Hence, the FFT is calculated at every dwell time (t s ) for N d times. 
4) Select components near DC according to 
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where K is the index of low pass filter cut off frequency (Bf) in the FFT, i.e. K = 
[N ■ Bf/F s \. This matrix reduction is in effect a low pass filter, that selects spectral 
feature of the primary signal and reduces noise power. The effective SNR is increased 
accordingly as long as the selected signal portion contains more power than the rest does 
in the same bandwidth. Algorithmic complexity is also reduced by a factor of N/ (2K + 1) 
due to the smaller matrix size in the covariance calculation. 
5) Calculate sample covariance of M as 



C = cov(M) 

= E fc [(M-l JVd M) T (M-l iVd M)] 



(6) 



where M = fi M = [/i , fj, 1 , ■ ■ ■ , //jv d -i], /V = m T 



i 



2K+1 



is the all one's vector with length N d and E fc [-] = ^ X^J']- Note that the covariance matrix 
C is symmetric. 



6 



6) Compute test statistic T = T x /T 2 , where 

^ N d -lN d -l N d -l N d -1 

Tl = aT E E c ™ = T * + aT E E c ™ < 7 > 

d r =0 m=0 r=0 u=t+1 

T 2 = W 52*r, (8) 

d r=0 

where c TU is the element of C at the r-th row and w-th column, which is the covariance 
of m T and m u . In other words, T 2 is the sample mean of the diagonal terms of the C 
matrix, which are the autocovariances of the spectrograms, and 7\ is the sum of all the 
elements of the covariance matrix. 

7) Compare T with decision threshold 7 as 

T ^ 7, S.t. 7 = arg SUp P F a{T) = PFA,req, (9) 

n r 

where Pfa(T) is the probability of false alarm with threshold T and PFA,req is the required 
false alarm probability given by the specification. If T exceeds 7, detection is declared 
(Hi ). Otherwise, the decision is made that no signal is present (Hq). 
The test statistics T may superficially look similar to [6], but there are several key differences 
that lead to SCS's better performance. First, J6l used correlation in the time domain while 
SCS uses the frequency domain. As shown in [|2) and llT6ll . power spectrum is one of the most 
significant features of the primary signal to be detected and is a frequency domain measure. 
Note that SCS detects unique spectral features that are not entirely flat, e.g. peaks or valleys, 
while [|6|| can be applied to any modulated signal without prior information. Second, SCS uses 
only part of the spectrum by low pass filtering in step 4), which results in reduced noise power, 
in other words, increased effective SNR. Downsampling in step 2) also helps by averaging out 
noise which has zero mean. Low pass filtering in step 2) is essential to prevent aliasing and 
suppresses noise power in stop band. Engaging this step can improve the detection performance 
of O considerably. Third, (6l uses non-negative (i.e. absolute) values of the covariance matrix 
but SCS uses both positive and negative values. 

III. Performance Analysis 

The performance of spectrum sensing can be measured by two probabilities: probability of 
detection P D and probability of false alarm Pfa, which are defined as 

P D = Pt(T> 1 \H 1 ), (10) 
P FA = Pr(T> 7 |H ). (11) 
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The primary goal of the proposed sensing algorithm is achieving high Pp and low Pfa at 
the lowest received signal power level, however these two probabilities trade off each other 
depending on the decision threshold 7. Therefore the required detection capability is determined 
by the application as a minimum P D and maximally allowed Pfa P a ir at me required SNR. 

In this section, detection performance of the SCS algorithm is analyzed according to the 
following steps. Note that the current signal model is one of the most difficult cases for the SCS 
detector, since it has only one peak over a wide flat spectrum. The analysis given in this section 
can be easily applied to other signal models with other spectral features, by suitably adjusting 
the value(s) of f c , K, Bf, t s , and N^. We will first set up the analytical model and derive the 
statistics of variables. Second, we will find a solution for Pfa- Since the noise is assumed to be 
white Gaussian and the Pfa analysis needs only the noise spectrum and its correlation, we can 
get the analytical results accurately. Then, the decision threshold 7 is obtained by setting the Pfa 
as required. Finally, the P D is analysed with the given threshold and received SNR. In order to 
find the theoretical solution for P D , the statistical distribution of T = Xi/T 2 is needed, however 
it is almost impossible to get in practice, especially when the primary signal can have arbitrary 
spectrum shape as allowed in the SCS algorithm. Therefore, we assume Gaussian distribution 
thanks to the central limit theorem as P F a case. Besides the Gaussian assumption, we simplifies 
some of the statistics for very low SNR regime that enables to find a closed form solution. A 
summary of the notations used in the analysis is given in Table HI The convention used in this 
paper is that the subscript p, s and w stands for pilot, the rest of the signal and noise respectively. 

A. Analytical Model and Statistics 

The partial matrix M captures spectrum on the band of interest, such as the spectrum near the 
pilot for the ATSC signal. Ideally, white noise would have a flat spectrum over the frequency 
band; however this cannot be true in a practical system using periodogram estimation. Hence the 
noise spectrum can be estimated by some fluctuations over the nominal noise power spectrum 
density (PSD) N w /2. The residual components, i.e. the fluctuations over the nominal spectrum, 
are weakly correlated due to the filtering, however the amount is small compared to the signal 
correlation. Since the diagonal elements of the covariance matrix C are the autocovariances, they 
are larger than the off-diagonal ones. When there is no signal components (Ho), C becomes 
almost a diagonal matrix and the resulting test statistic T gets close to 1 . When there is a primary 
signal whose spectrum is not flat, the residual spectrum is highly correlated and T increases. 



When there is no signal, the spectrum is estimated as 
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The signal and noise components are assumed to be uncorrelated. 
If there is no signal, the statistics of m T (k) are given as in ||2TT| 

E[m T (k)} = N /2 = N w /2 4 ^ 
var(m T (A;)) = (N w /2) 2 , 



(14) 



(15) 
(16) 



where iV is the one-sided PSD of the noise signal w(t), which has full span of the bandwidth 
B. Note that fx T E[m r (fc)] for non-white signal. Since w(n) is the filtered and downsampled 
version of w(t), the effective noise power is given as 

a 2 w = a 2 Q B f /B = N w B f , (17) 

where a 2 Q is the noise power in received signal x(t) = w(t). A perfect low pass filter with cutoff 
frequency Bf is assumed here. 

In this paper, we assume the primary signal has a pilot tone, which has f p (0 < f p < 1) of the 
total signal power Pg. With these settings, we can obtain the mean PSD of the received signal 
as follows. 

Lemma 1: The in-band average spectrum of the signal and noise at the r-th sensing window, 
which is the mean of the r-th column of the matrix M, is 



/'t 



N w + 



<Ut) N s {t) 



2K 



(18) 



Proof: See Appendix \M 
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In the low SNR regime where N w ^> N s , \x T can be further simplified as 

SJr) 

Vr*N w + J^L, (19) 

i.e. the in-band power is mostly concentrated on the pilot frequency and it can be comparable 
to the noise power even in the extremely low SNR where the rest of the signal power is almost 
dominated by the noise components. 

The elements of the covariance matrix C can be obtained as 



-"UT 

K 



E k [(m T - /i T l 2 K+i) T (m u - /i„l 2 K+i)] (20) 
1 K 

2K ( mr ~ ^t) (m* (k) - (21) 

k=-K 

= E fc [m T r m u ] - fi T fi u ( ^ K ^ - j • (22) 

Define normalized correlations for signal power spectrum and noise as 

, _ E[b p (r)b p (u)] E[N s (t)N s (u)] _ , 
Up - „ > N2 - a s , ^) 

u p s 

E fc [(d;)^d^] = ^J2 k S w (r,k)S w (u,k) 
VWd£)vax fc (d£) N* 
where I = t - u and d T w = [5 w (r, —K), ■ ■■ } b w (r, K)] t . 

Then the statistics of c ru cun be obtained using the following Lemma [2l which specifies the 
covariance matrices. 

Lemma 2: When the primary signal is present, the statistics of c TU can be derived as 

a 1 1 
E[c TU ] = a l w Nl + ^5 2 p + N s (N w + (N s + N w )a l s )+0(—), (25) 

+0(^j), (26) 

where I = r — u and bin) is a Kronecker delta function, which is bin) = { o If n^o • 

Proof: See Appendix [Bj ■ 
Since K ^> 1, the last terms of (1251) and (1261) can be omitted for simplification. As in (fT9l) . 
the expectation of be further simplified in the low SNR regime as 

E[c TU ] « aJX + |p p 2 . (27) 
Note that the no signal case "Ho can be found by letting N s = b p = 0. 
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The Lemma [2] clearly shows that the covariance matrix depends on the correlation of the 
received signal and its relative power to the noise power, which is the SNR. The statistics of 
the test are given in the following Theorem [TJ 

Theorem 1: The statstics of the test crieterion in the low SNR regime are 

E[Ti(JV d )] « N l + ^ + 2 -3± N l + §^5l (28) 



2K N d w KN t 



E[T 2 (N d )} » A^ + JL, (29) 

if /y2 ^ 

var(T 2 (iV d )) » |^(A^ + <y 2 + -f {N 2 W K + (N w + 5 p ) 2 )j , (30) 

where the accumulated correlation of the noise and signal are given as 

N d -2 

A N d = £ - 1 - l)oL, (3D 

1=0 

N d -2 

A N d = Y,(Nd-l-lH. (32) 

1=0 

Again the no signal case is obtained by letting all the signal terms to be zero and the approxi- 
mations become exact. 

Proof: Equations (1281) through (1301) can be easily derived from the results of Lemma [2] and 
d27l) as 



E[T 2 ] = E[c TT ], (33) 

iV d -2 N d -1 

E[Tx] = E[T 2 ] + — ^ ^ E[c r J (34) 

d T =0 u=r+l 

var[T 2 ] = -^-var(c rr ). (35) 

U 

The approximations are made here based on the simplified statistics of the covariance matrix 
for the low SNR regime, and will be verified through simulation in Section [IV] The results of 
the Theorem Q] shows that T 2 is composed of autocovariances of the signal and noise term and 
T\ is composed of accumulated covariances. Therefore, we can expect the relative amount of 
these values will determine the detection performance, which will be derived in Section IIII-Bl 



B. Probability of Detection 

In this section, the performance of the SCS algorithm is analyzed in two steps. The decision 
threshold is found first and then the probability of detection is derived as in J6]|. By the central 
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limit theorem, both T\ and T 2 approach Gaussian distribution when the number of samples 
(N d x N) becomes large. The decision threshold and the corresponding probability of detection 
are given in the following Theorem |2] and [3j 

Theorem 2: The decision threshold 7 is determined as 



7 



2v ^(l + ^) 



Q-^l-P FA ^ q ) + 2^/KN d 
where PFA,req is the required probability of false alarm. 

Proof: The probability of false alarm P FA can be calculated as 
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where 
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Therefore, 7 is easily derived from © and d38l) . ■ 

If we assume the noise correlation is independent of the lag, i.e. a l w = a w ,Vl, then A^ d 
becomes proportional to iVj for large values of K ■ N d product and the threshold 7 has a 
tendency to become proportional to N d . Note that the decision threshold is independent of the 
noise power and SNR, which explains the robustness of SCS to the noise power uncertainty. 

Theorem 3: The probability of detection is obtained as 
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where T p is the effective SNR, which is 
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Proof: The probability of detection can be derived in a similar way as Theorem [2l The 
variance of T 2 can be further approximated for the low SNR case as 



var(T 2 (iV d )) ^ i- (| (N w + 5 P ) + ^ (n w K + | 
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Then, the probability of detection is derived the same way as the proof of Theorem |2] as 
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/ 



Given SNR and sampling rate F s , the effective SNR is determined by the number of FFT iV 
in a dwell, since the signal bandwidth B and pilot power fraction f p are fixed for the ATSC 
signal. Theorem [3] indicates that the detection probability Pp can be increased by increasing the 
number of samples either by N or N^. Note that K is also increased when N is increased. Also 
note that the decision threshold 7 is increased along with Nd- Hence, we expect that the length 
of sensing window or dwell t s affects the performance more than the number of dwells N d . 



IV. Simulation Results and Comparisons 

In this section, we present simulation results of the proposed SCS algorithm that show its 
practicality and superior detection performance versus existing solutions in terms of sensitivity 
and sensing time. The simulation has been conducted in compliance with the IEEE 802.22 
standards using actual DTV signals captured in the US, for fair comparison and more importantly, 
to prove SCS's applicability to real spectrum sensing problems. Our simulation codes are publicly 
available in lfl9l . The analytical results driven in Section |In] are also verified in this section. 
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A. Spectrum Sensing Requirements 

In any CR based system, emphasis is on the sensing of potential primary services in the given 
frequency band. Each CR system has its own requirements on the spectrum sensing depending 
on various aspects of the primary services, propagation characteristics and application specific 
features. Spectrum sensing requirements of the IEEE 802.22 system are reviewed in this section 
and used in the simulation. 

The sensitivity is defined as the minimum received power at the secondary system that can 
be detected with the given probability of detection (Pd) and probability of false alarm (Pfa) 
pair. Note that the sensitivity is not equal to the SNR wa n in [[T4|. The primary services in the 
IEEE 802.22 systems are digital TV, analog TV, and wireless microphone, operated in TV Bands 
(54 - 862 MHz) Q, [|22]|. Note that the format and bandwidth of the TV signal depends on 
the regional regulation. For example, digital TV signals can be Advanced Television Systems 
Committee (ATSC: North America), Digital Video Broadcasting - Terrestrial (DVB-T: Europe) 
or Integrated Services Digital Broadcasting (ISDB: Japan) 0, [|9]|. hi this paper, only ATSC 
A/74 DTV signal Il23l is considered for detection. The required sensitivity is 90% of Pd and 
10% of Pfa at -H6 dBm, which is equivalent to -21 dB of SNR. Figure [T] shows a sample 
ATSC spectrum that will be detected using the proposed algorithm. As explained in Section |inj 
it has a high powered pilot on the lower edge and is correlated near the pilot. 

B. Simulation Procedure 

The simulations were performed according to the procedures in (|24| and 11251 with the 
proposed SCS algorithm. Two previous sensing algorithms are also tested for comparison and 
the performance analysis given in Section [III] is also verified. The actual DTV signals captured 
by the ATSC in Washington DC and New York City are used for simulation, which includes 
all relevant impairments of wireless channels such as fading, scattering, shadowing and local 
oscillator mismatch. The 12 selected signals recommended by the 802.22 committee (|24| are 
used as in several other studies 0, EH, lfT5l . 

Due to the random characteristics of the wireless channel, detector performance varies signif- 
icantly. Figure [2] shows the probability of missed detection (Pmd = 1 — Pd) f° r five selected 
signals that illustrates the broad range of sensing performance along with the average of the 12 
signals. Note that the sensitivity difference can be as large as 12 dB with Pfa = 0.1. The 21.52 
MHz sampled ATSC captured signal is downsampled to F s = 2.152 MHz. The FFT point N is 
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determined by 

N = T s.t. n = Llog 2 (F s /A/)J = Llog 2 (F s • t s )\ , (47) 

where A/ is the frequency resolution of the FFT and t s is one sensing time. For Fig. |2l t s = 1 
ms and N = 2048 are used. The number of dwells N d is 30 and LPF bandwidth Bf is set 
to 20 kHz to provide a buffer for frequency deviations. The sensing requirements of 802.22 
are achieved with these parameters. Downsampling frequency F s and bandwidth Bf are fixed 
throughout the rest of the paper. 

C. SCS Sensing Results and Comparisons 

Before we evaluate the probability of detection, the analytical value of decision threshold is 
verified in Table [H] with simulation. The decision threshold 7 is roughly proportional to A^ as 
expected. It is slightly reduced as t s increases, but the effect is negligible compared to that of 
Nd- Also note that the analytical value fits better as the number of samples increases as expected. 

The detection performance of the proposed SCS detector is compared with the FFT based 
pilot location detector 0, covariance absolute value (CAV) detection Q, and CAV with low 
pass filtering and downsampling (step 2 of the SCS algorithm) in Fig. [3] on detecting ideal ATSC 
signal with additive white Gaussian noise (AWGN). Sensing times for SCS and the pilot location 
detector are set to 30 ms (1 ms x 30 and 5 ms x 6 each) and CAV used iV s = 5x 10 5 samples 
(23.2ms) with a smoothing factor L = 14. The decision thresholds for SCS and CAV are set 
to make Pfa = 0.1. The P FA of the pilot location detector was 0.02. The sensitivities of SCS 
and pilot location detector are about 7 dB better than that of the CAV, while CAV with step 2 
outperformed SCS about 1.3 dB. Note that the step 2 not only reduces the complexity but also 
significantly improves the sensitivity of the CAV detector due to the increased effective SNR. 
This step is also used in J21- 

Though Fig. [3] shows that SCS, pilot detector and improved CAV meet the sensitivity re- 
quirement of the IEEE 802.22 standard, it does not guarantee their performance in the real 
environment. Since detector performance is highly dependent on the received signal quality that 
is degraded by the unpredictable wireless channel, it is very important to test the detector under 
practical circumstances. Hence, all the simulation results from now on are conducted using the 
actual signal sets explained in Section IIV-B[ 

Figure 0] shows the superior sensitivity of SCS based on the number of dwells A^. All the 
parameters for each detector are set as used in Fig. [3l The sensitivities are averaged over all 12 
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signals for each detection method^ The SCS algorithm shows similar performance with pilot 
location sensing when N d = 6, which is l/5th of the sensing time. Fast sensing has several 
advantages. It can reduce interference to the primary service by enabling fast channel changes 
ifTTI . Furthermore, the secondary system can use the channel more efficiently since it cannot 
transmit during the sensing time E6l . When the same sensing time (30 ms) is used, it shows 
about 3 dB better sensitivity. Figure |4] also shows that SCS is very effective detector in real 
environment. Though improved CAV (with low pass filtering and downsampling) has better 
sensitivity for ideal ATSC signal than SCS, the SCS outperforms it when detecting actual signals. 

Given the sampling rate, the frequency resolution or equivalently number of FFT size iV is 
determined by the duration of the sensing window t s as in (PTTl) . For example, 0.1 ms sensing 
window allows a 128-point FFT, while 1 ms allows a 2048-point FFT. In Fig. [5l we show that 
fine frequency resolution (a longer sensing window) detects better when the same dwells are 
used, which also coincides with the analytical results. The sensitivity of the SCS algorithm is 
improved by roughly 2 dB per doubled window size. 

Combining the above results, the relationship between the total sensing time T s = t s x Nd 
and detection performance of SCS is investigated in Fig. [6l The simulation results show that 
t s affects performance more than N d , as expected from Theorem [3l Note that 12 ms sensing 
time with t s = 2 ms has sensitivity of -19.9 dB, while 15 ms with t s = 0.5 ms has -19.3 dB. 
It clearly shows that there is a tradeoff between performance and complexity. Increased sensing 
time allows an increased FFT size. One way to reduce complexity is to decrease the down 
sampling frequency F s and to use coarse frequency resolution. However, since a 2048-point 
FFT is included in every 802.22 device, it can be reused in spectrum sensing. Thus, using a 1 
ms sensing window will be a reasonable choice. Sensitivity of -22.9 dB is achieved with 2 ms 
x 30 = 60 ms sensing time. 

The noise power has thus far been assumed to be known at the receiver. However, perfect 
estimation of noise power is impossible in practice and it is well known that even a fraction 
of dB noise uncertainty can degrade detector performance severely 031. Figure [7] shows the 
robustness of SCS to the noise uncertainty problem. The noise uncertainty factor p means that 
the noise power has a uniform distribution of [— pdB,pdB]. SCS shows almost no performance 
degradation with 2dB of noise uncertainty regardless of sensing time. Note that the nominal 



2 We corrected the average calculation in (2): it was done with 9 good signals. 
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noise uncertainty level for the IEEE 802.22 system is 1 dB [22J. 

The robustness of SCS mainly comes from the fact that it exploits the statistical independence 
of the signal and noise components, especially that the noise is uncorrected. Thus, uncertainty 
in the noise only reduces the possible correlation and does not strongly affect signal detection. 



In this paper, we have proposed a spectrum sensing algorithm using the covariance of the partial 
spectrum of the received signal. The proposed SCS algorithm can be used to detect arbitrary 
signals with a proper selection of parameters, depending on the features of the primary signal. 
The decision threshold in particular should be carefully chosen since there is a fundamental 
tradeoff between the probability of false alarm and the probability of detection. We derived an 
appropriate decision threshold and the detection time and sensitivity of SCS mathematically. The 
theoretical results were verified through extensive simulations conducted using actual digital TV 
signals captured in the US, and our code is openly posted online [19]. The detection performance 
of SCS is compared with the FFT based pilot location detector and the CAV detector. The SCS 
detector achieves the same detection performance as the FFT based pilot location detector in 
l/5th of the sensing time and 3 dB better sensitivity in the same sensing time of 30 ms. It has 
also been shown that SCS is highly robust against noise uncertainty. The results of this paper 
suggest that SCS is an extremely effective sensing algorithm and should be seriously considered 
by the 802.22 standards body, and future industry efforts at cognitive radio. 



The spectrum can be categorized into three sections: noise only (lower frequency, k ), pilot 
tone (fci), and signal data part (high frequency, k 2 ). Without loss of generality, we can set K t 



V. Conclusion 



Appendix A 



Proof of Lemma [j] 



at DC, i.e. K 



{-K,-K+l r .- 



1}, Ki = {0}, and K 2 = {1, • • • 



K}. 



Then, 



E[S w {t, k)] 



E[S s (t, k)] = 0, 



(48) 



E[m T (k)} 



N w + S p l Kl (k) + N s l K M 
(N w + 5 p l Kl (k) + N s l K2 (k)) 2 



(49) 



var (m T (k)) 



(50) 
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where 5 P = E[5 p (t)) is the pilot power spectrum, N s = E[iV s (r)] is the nominal PSD of the 
primary signal and 1a (k) is an indicator function. The in-band average spectrum is given as 

= ml (N w + 5 w (r,k))+ (N w + 5 w (r,k) + 6 p (r)) + 
( N ™ + S ™( T > k) + N s (t) + 8 s (t, k)) 

k£K 2 

where 5 P = ^ = jfihPs, N s = (1 ~^ )Pg and Af = l/t s = F s /N is the frequency resolution. 

Appendix B 
Proof of Lemma [2] 

The expected value of be easily shown using ([HI), CGI and <@8]) - (ED- When there 

is signal, the first term of (|22j) can be analyzed in three frequency regions as 

1) keK : 

E[m T T m u ] = E[{N w {r) + 5 w {r,k)){N w {u) + 5 w {u,k))} 

= N 2 + N 2 w a { r u) = N 2 W (1 + oL) , (52) 

2) keKi. 

E[m r T m u ] = E[(N w (r) + 5 w (T,k)+5 p (T))(N w (u)+5 w (u,k)+5 p (u))} 

= N 2 (l + a l w )+2N w 5 p + a l p 5 2 p , (53) 

3) keK 2 : 

E[m T T m u ] = E[ \N w (t) + 5 w (t, k) + N s (t) + S s (r, k) + ^Re{S T (k)W:(k)}^j ■ 

N w {u) + 6 w {u, k) + N s (u) + 5 s (u, k) + ^Re{S u {k)W:{k)}\ } 

= E[N w (t)N w (u)} + E[6 w (t, k)5 w {u, k)] + 2N W N S + E[N s (r)N s (u)] 

+E[5 s (r)8 s (u)] (54) 

= N 2 (l + a l w )+2N w N s + 2N s (N s + N w )a l s , (55) 

E[Re{S T (k)W*(k)}] = E[Re{S T (k)}Re{W T (k)}] + E[Im{S T (k)}Im{W T (k)}] = 0, 
E[5 2 s (t, k)] = var (m T (k)) - E[5 2 w (t, k)] = N 2 + 2N W N S . 
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The expected value of the second term of (1221) can be evaluated straight forwardly 



E[/i T /iJ 



E 



2K 



N w + 



2K 



, N w 6p 



+ 



a 



p sl 



K ' AK 2 p ' 

Now the expectation of the sample covariance is given as 



(56) 



E[c T 



IK + 1 r, / , x 5n 

-N 2 W 1 + at,) + 



2K 
2K+1 



2K + 



j(2N w + a%) + N s (N w + (N 8 + N w ) a l s ) 



2K 



2 , Nw°p 

w K 



+ 



a 



P X 2 



[K 2 S P 



J 



ulK + ^5 2 p + N s (N w + (N a + N w )a l s ) + O(-) 



1 

K 2 ' 



Next, we will derive the variance of c TU for the signal case. Let 

V T m (k)=m T (k)-fi T . 
Then the variance of be obtained as 

' ' ~ l ' ^var(V7(*)V7(*)), 



var c 7 



\ fe / k 



(57) 



(58) 



(59) 



where V T m (A;) and V^ TO (Jfe) are i.i.d and V™{k)V™{k) and y™(j)K m (j) are uncorrelated if Jfe ^ j. 
Since the variance of the product of two random variables X and Y is given by 

var(XF) = (E[X}) 2 a 2 x + (E[y]) 2 tr£ + 2E[X]E[y]cov(X, y) + o\ + erf- + (cov(X, Y)f , (60) 

where cr^ and a\ are the variances of X and y respectively. We need to find the statistics of 
V™(k) and V™{k) first. It can be shown that 



E[v T m (k)} = nv™(k)} 



2K' 



'pi 



k e K , 
k e K u 



N s -^, keK 2 , 



var(y T m (A;)) = var(y™(A;)) = var(m r (fc)) = (N w + 5 p l Kl (k) + N s l Ki (k)f 



(61) 



(62) 



cov (V™{k), V™{k)) = var(V™(k))5(T - u). 



(63) 
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Now we can find the variances in (l59l) . Define Cf = v&r((V™(k)V™(k)) , k e kj, then 



E 



2K 



1 2 



N„. + E 



£ p (w). 



2iT 



TV 



-2E 



E 



S p (u) 
' 2K 



N 2 J{r - U ) + N 2 W + N 2 W + Ni5{r - u) 



= K (JL (1 + a l p 6(r -u))+2 + N 2 J(r - «)) + 0{±), 
G\ = 2(N w + S 2 p ) 2 (S 2 p (l + a l p S(T-u)) + (N w + 5 p ) 2 5(r-u))+0(^) 



2N 2 



5, 



p - (1 + a l p 5{r - «)) + (iV w + iy s ) 2 5(r - «) ) + 0(— ). 



AK 2 

Thus, the variance of c tu C3.n be obtained from (l59l) as 

, , 1 



var c T 



(64) 
(65) 
(66) 

(67) 



when there is a detected signal. The no signal case H is obtained by letting N s = 0. 
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TABLE I 

Notations used in Analysis 



Variable 


Description 


Unit 


s T (n), w T (n) 


signal and noise samples on the r-th dwell, s r (n) = s(n + tN) 


V 


m T (k) 


k-th component of the spectrogram vector on the r-th dwell m T 


W/Hz 


N s (t), N w {t) 


average signal and noise spectrum computed in the signal band at the r-th dwell 


W/Hz 


S s (t, k), 5 w (t, k) 


residual signal and noise components at the fc-th bin and the r-th dwell 


W/Hz 


S P (r) 


power of the pilot tone at the r-th dwell 


W/Hz 


S T (k), Wr(k) 


fourier transforms of s T (n) and w T (n) 






in-band average spectrum at the r-th dwell 


W/Hz 


N = N W 


nominal one-sided PSD of the noise signal 


W/Hz 


N s 


nominal signal PSD 


W/Hz 


S P 


nominal pilot power 


W/Hz 


Ctu 


sample covariance of m T and m u 




„l t 

a p , a s 


nomalized correlations for pilot and signal each 






nomalized correlation for noise term 




AP A w 
A N d ' A N d 


accumulated correlations of pilot and noise 





TABLE II 

Decision Thresholds for ATSC signal detection with P F a = 0.1 



Analytical Value <f36t 


N d 


t s (ms) 


0.1 


0.5 


1.0 


2.0 


6 


3.39 


2.83 


2.73 


2.70 


12 


5.92 


4.89 


4.70 


4.64 


30 


13.16 


10.96 


10.60 


10.36 


Simulated Value 


N d 


t s (ms) 


0.1 


0.5 


1.0 


2.0 


6 


3.59 


3.07 


2.92 


2.75 


12 


6.51 


5.27 


4.29 


4.79 


30 


14.17 


11.27 


10.95 


10.62 
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Fig. 1. Spectrogram of a captured ATSC signal around pilot tone at 2.69 MHz. 




Fig. 2. Probability of missed detection of Spectral Covariance Sensing (SCS) for the 5 ATSC signals selected from the 12 
signals set, t s = 1ms, = 30. 
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Fig. 3. Probability of missed detection of Spectral Covariance Sensing (SCS), FFT pilot location detector |2|, Covariance 
Absolute Value (CAV) [6], and CAV with low pass filter and downsampling (step 2 of SCS algorithm) on the ideal ATSC signal. 




-30 -28 -26 -24 -22 -20 -18 -16 -14 -12 -10 

SNR(dB) 



Fig. 4. Effects of Nd in Spectral Covariance Sensing (t s = 1ms) compared with the FFT pilot location detector (2) and 
covariance absolute value (CAV) detector |6|. 
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SNR(dB) 

Fig. 5. Effects of t s in Spectral Covariance Sensing, N d = 12. 




SNR(dB) 



Fig. 6. Probability of missed detection of Spectral Covariance Sensing in various sensing times (t s x Nd). 
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Fig. 7. Comparison of the probability of missed detection of Spectral Covariance Sensing with 2 dB noise uncertainty and 
dB (without uncertainty). 



